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Abstract 

The deformation energy of the even-even nuclei of the Cerium isotopic chain is investigated by means of the Macroscopic- 
Microscopic method with a semiclassical shell correction. We consider axially symmetric shapes. Binding energy and two 
neutron separation energy are also evaluated. For the sake of clarity several important details of the calculations are also given. 
It turns out that all these nuclei have prolate equilibrium shape. The regions of maximum deformation are obtained around 
N = 64 and N = 102. There is no critical-point of quantum phase transition in this isotopic chain. 
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I. INTRODUCTION 


Nowadays it is well established that the majority of nuclei possess a nonzero intrinsic electric quadrupole moment 
(IEQM). This feature means that the charge distribution inside the nucleus deviates from the spherical symmetry. 
In other words, apart from very few nuclei, the surface of the nucleus is generally not spherical in its ground state. 
The intrinsic quadrupole electric moments (or equivalently the nuclear deformation) can be deduced from two types 
of measurements: 

• The reduced electric quadrupole transition probability, B(E2) jlj 

• The static electric quadrupole moments of ground and excited states, Q 0 

It turns out that in a number of cases, the two methods of measurement do not systematically lead to the same values. 
Important discrepancies occur for several nuclei. This is essentially due to the fact that not only different experimental 
techniques are used but above all, because different models can be implemented to deduce the nuclear deformation 
for the both cases. In Ref.Q it is stated that deformations deduced from B(E2) have a ’’more general character”. 
In other words, “i?(£ , 2)-type” data reflect not only static nuclear deformation (permanent deviation of the nuclear 
shape from sphericity), but also dynamic deformation. Furthermore, B(E 2) measurements are model independent 
and thus are generally more reliable. This is corroborated by the fact that the only systematic compilation in which 
the deformation of the ground state is given explicitly is based on B(E2; 0 + —>• 2 + ) and has been published in Ref.[l[. 
In the present work, experimental values refer to these ones. 

Theoretical approaches to the deformation energy can be divided into two categories; Dynamic calculations to 
find the shape of the ground state (or even of excited states) and static calculations by determining the absolute 
minimum (ground state) or multiple minima (shape isomers) in the potential energy surface (PES) for a given nucleus. 
Thus, on the one hand, we have the so-called collective models, which themselves are subdivided into two groups: 
The ’’Geometric Collective Model” also called ’’the Collective Bohr Hamiltonian” (CBH and its variants) and the 
’’Algebraic Model”, well known under the name of the ’’Interacting Boson Model” (IBM and its variants) 0. On the 
other hand, ’’particle models” consider the nucleus as a collection of interacting nucleons (fermions). In practice, the 
classical JV-body problem can be approximately solved by the usual approximation of the mean field with eventually 
residual interactions. In this respect, the “best” mean field is deduced after applying a variational principle in the 
Hartree-Fock-Bogoliubov method (HFB). In this model, the determination of the potential energy surface (PES) of 
the nucleus amounts to perform constrained Hartree-Fock-Bogoliubov (CHFB) calculations Q. We will not address 
very complicated methods “beyond the mean field” such as the Quasiparticle Random Phase Approximation (QRPA) 
or the Generator-Coordinate-Method (GCM) methods which are unsuitable in practice for large scale calculations. 
Because of CHFB calculations are time consumers, especially in large studies, Microscopic-Macroscopic method (Mic- 
Mac) constitutes a good alternative which, is up to now, implemented 0. In the present work, we use an improved 
variant of this method. The word “improved” means that we use semi-classical method to avoid the well-known 
drawbacks (spurious dependence on two mathematical parameters) of the standard Strutinsky shell correction (see 
text below). 

The present study is devoted to the deformation energy, equilibrium nuclear shapes and binding energy of the 
ground state of the even-even cerium isotopes. There are many reasons to this choice. One of them is to re-test our 
previous calculations. In effect, similar calculations have been already performed by us in the xenon, barium, and 
cerium region Ref [7|. However because the phenomenological mean potential varies smoothly with N and Z , we have 
made, in the past, a rough approximation by choosing the same set of parameters for the phenomenological mean 
potential, for the all treated nuclei. Originally, this approximation was done only for simplifying the calculations. 
Here, contrarily to that study, each nucleus has its “own” mean potential with a specific set of parameters. In this way 
it is possible to evaluate in a rigorous way the uncertainty introduced in the previous calculations. Apart from this 
remark, there are several main other reasons which could justify this choice: (i) First, it should be interesting to see 
how the deformation energy and binding energy vary with the neutron number (N) for this isotopic chain, (ii) Second, 
the present study extends the previous calculations to all cerium isotopes up to the drip lines (34 versus 13 nuclei), 
(iii) Third, we also will attempt to deduce, from potential energy surface (PES) curves, the shape transition from 
spherical to axially deformed nuclei, looking for the so-called A'(5) critical-point between U(5) and SU(3) symmetry 
limits of the IBM 1-0. 

It is worth to recall briefly some information deduced from the literature for the cerium isotopes. In the past, a 
number of experimental as well as theoretical studies have been done for the cerium isotopes. Among the numerous 
studies, we only cite some of them: In 2005 Smith et al 0 have studied excited states of 122 Ce up to spin 14 h 
deducing a probable quadrupole deformation of about /? ~ 0.35. The deformed nucleus 130 Ce has been studied in 
1985, using the techniques of in-beam gamma -ray spectroscopy El- The corresponding data have been interpreted in 
terms of the cranking model by assuming a prolate deformation with £2 ~ 0.25 (/3 ss 0.27). High-spin states in 132 Ce 
have been also studied by A.J. Kirwan et al. They found a superdeformed band with deformation /3 « 0.4 much 
more larger that the ground state deformation (/? « 0.2). E. Michelakakis et al 0] by evaluating 7 —ray transitions in 
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142 Ce and 144 Ce conclude that in cerium isotopes (near the beta-stable line) the onset of nuclear deformation occur 
between N = 86 and N = 88. ’’Pure” theoretical calculations have been performed in Ref. 0 and 0 with projected 
shell model (PSM) and Hartree-Bogoliubov ansatz in the valence space respectively for 122 Ce and 124 ~ 132 Ce for low 
lying yrast spectra. Good values of energy levels and reduced transition probabilities B(E 2,0 + —> 2 + ) have been 
obtained respectively in these two papers. Other approaches for the rich-neutron cerium isotopes have been made 
in Ref. 0- A study of the shape transition from spherical to axially deformed nuclei in the even Ce isotopes has 
been done in Ref |17j using the nucleon-pair approximation of the shell model. The result of a such study is that the 
transition has been found too rapid. Relativistic Hartree-Fock-Bogoliubov theory has been used to predict ordinary 
halo for 186 Ce, 188 Ce, 190 Ce, and giant halo for 192 Ce, 194 Ce, 196 Ce, 198 Ce near the neutron drip line. Systematic 
studies about nuclear deformations and masses of the ground state can be found in Ref. 0-0 with respectively, the 
Finite-Range Droplet-Model (FRDM) , Hartree-Fock-Bogoliubov (HFB), HFB+5-dimensional collective quadupole 
Hamiltonian and Relativistic Mean Field (RMF) models. 


II. THE MACROSCOPIC-MICROSCOPIC METHOD 
A. Liquid drop model and microscopic corrections 

This method combines the so-called semi-empirical mass formula (or liquid drop model) with shell and pairing 
corrections deduced from microscopic model. Thus the binding energy is given as a function of nucleon numbers and 
deformation parameter (referred to as /3) by mean of the usual symbols: 


B(A, Z, P) = E ldm {P) - SB micro ((3) (1) 

SB micro contains the shell and pairing correction (see text below). The minus sign before 6B m i cro is consistent with 
the convention that the binding energy is defined as positive here. For the liquid drop model we take the old version 
of Myers and Swiatecki |28) (because of its simplicity compared to more recent formulae). Here, there is no need to 
look for very high accuracy in binding energy, because this is not the purpose of the present work. 

Eldm(P) = C V A - C s A 2 / 3 BsU 3) - C c Z 2 A~ 1 / 3 B c {P) + ea pair A- 4 / 2 + C d Z 2 A~ 1 (2) 

In Eq.fl5]), we have the usual contributions of volume, surface and coulomb energies. The different constants of Myers 
and Swiatecki are given in appendix A. The shape dependence (/3) of the surface and coulomb energies are contained 
in Bs(/3) and Bc(/3). They are normalized to the unity for a spherical nuclear surface. The latter is symbolized by 
(3 = 0. The two last terms in Eq.@ are respectively due to the smooth part of the pairing energy and the correction 
of the Coulomb energy to account for the diffuseness of the nucleus surface. The different constants will be fixed later. 
The potential energy surface (PES without zero point energy correction) is defined as follows: 


Epes{P) = E ldm { 0) - B(A, Z, 0) = A E ldm (P) + SB micro {/3) (3) 

in which 

A E ldm {/3) = E ldm ( 0) - E ldm {/ 3) = CsA 2 A [Bs((3) - -Bs(O)] + CcZ 2 A 1 ^ 3 [Bc(P) — -Bc(O)] (4) 

Constants Cy and Cs are expressed by means of three other constants ay, as, and n. For spherical shape, as said 
before, the normalization is expressed by: Bs( 0) = Bc{ 0) = 1. As it can be easily seen, the potential energy surface is 
related only to two macroscopic constants Cs (which depends actually on as and k) and Cc- To calculate microscopic 
shell an pairing corrections contained in SB m i cro , we have to proceed in two steps. The first consists in solving the 
Schrodinger equation and the second in deducing the shell and pairing corrections in an appropriate way, as explained 
in the following. 


B. Microscopic model 

We briefly present the microscopic model which is based on the Schrodinger equation of the deformed independent 
particle model: 

H(f3) I *,;(/?)} = £i (/3) I v^CS)) (5) 
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where | \Eq) and £j are respectively the eigenfunctions and the associated eigenvalues of nucleons. Hamiltonian H 
contains four contributions which are: (i) kinetic energy, (ii) central deformed mean field, (iii) spin-orbit and (iv) 
Coulomb interactions. We perform analogous calculations as in Nilsson model but our deformed mean potential 
is of Woods-Saxon type and therefore is ’’more realistic”. Although calculations are not self consistent, they are 
microscopic. It is to be noted that our Schrodinger equation has a form which is very close to the one of the Skyrrne- 
Hartree-Fock method. 

Eq.© is solved by our FORTRAN program described in details in Ref { 22 J and improved in two successive versions 
[H and 13. 


C. Microscopic corrections 

Microscopic corrections are defined as the sum of shell and pairing corrections which themselves are calculated 
separately for each kind of nucleons. 

SB m i cro (/3) = SE shell (N, /3) + 6E shell (Z, f3 ) + 5P pa i r i n g(N : /3) + 5P pa i r i ng (Z , /3) (6) 

In this formula the shell correction is defined by the usual Strutinsky prescription, i.e., as the difference between the 
sum of the single particle energies (which contains the shell effects) and an averaged (or smoothed) sum (which is free 
from shell effects) 


NorZ _ 

5E she u(N or Z) = £s<(ffl-I>OT (7) 

2 = 1 2=1 

Energies £j(/3) are deduced from Eq.©. In our procedure, the second sum is found by means of a semi-classical way 
instead a Strutinsky smoothing procedure, see Ref. [ 27 } • This avoids the well-known weakness of the standard shell 
correction method, namely, the dependence on two unphysical parameters which are the ’’smoothing” parameter and 
the order of the curvature correction. Moreover, it has been clearly shown that Strutinsky level density method is 
only an approximation of that of the semi-classical theory [26]. 

The ’’pure” pairing correlation energy is defined by: 


00 N/ 2 orZ /2 2 

P (P) = J2 2£i ^ V i ~ 2£i 0 9 ) ~ ~Q 

2—1 2=1 

where v'~ , A and A are the usual occupation probabilities, gap and Fermi energy of the BCS approximation (the factor 
”2” is simply due to the Kramers degeneracy). Since the smooth part of pairing correlations is already contained in 
the liquid drop model, we have to add only the one due to the shell oscillations of the level density. This contribution 
is defined by means of a similar formula to Eq. © 

SPpairingiN Or Z, (3) = P(/ 3) - P{0) (8) 

-- - 2 

where the averaged pairing is defined as P{/3) = (l/2)g sem i c i,(X)A . We use a simple BCS method to account for 
pairing correlations. To calculate Eq © and © we follow the method detailed in Ref.[27j with its FORTRAN code. 
The treatment of the pairing has also been explained in Ref. Q and references quoted therein. 


D. Numerical constants and prescriptions 

1. Constants of the Microscopic model 

For each kind of particles the mean central and the mean spin-orbit field are written as [13 : 


V(f3) 


_ Vo _ 

1 + exp (R v L v (P)/a 0 ) 




_ Vo _ 

1 + exp(RsoLso((3)/ a o) 


(9) 
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where Ly(/ 3) and Lso{P) contains the information on the deformation. In fact, these functions contain 9 constants: 
ho neut: ho prot^Jj^Vneuti Rvprot, Rso—neut, Rso—prot, c,(), A neuti A prot• These Quantities are taken from the universal 
parameters [291 ] (see appendix B) which is an optimized set. 

The Coulomb mean field is approximated by a uniform charge distribution inside a deformed surface. The volume 
conservation is therefore Vol = ( 4 / 3 ) 7 ri?^with the simple assumption R c h = Rvprot- 


2. Constants of the liquid drop model 

As already stated, we have chosen the parameters of Myers and Swiatecki (see table QJ because this set contains a 
reduced number of parameters with respect to more modern formulae. 

All the constants are needed in the binding energy whereas only as, Cc, k play a role in the potential energy surface. 


3. Nuclear mass excesses 

Nuclear masses are deduced as mass excesses: 

M excess (A, Z) = ZM h + {A- Z)M n - B(A, Z) 

where Mh = 7.289034MeU is the hydrogen mass excess and M n = 8.071431MeU the neutron mass excess. This 
makes comparisons with experimental values easiest. 


III. RESULTS 

In our previous paper Q calculations for isotopes 116 ~ 130 Ce showed that the equilibrium deformations (/3 ~ 0.25 — 
0.30) have always been obtained for symmetric prolate shapes (7 = 0°). Results obtained in Ref.[32j with a similar 
approach for the nuclei 116 ~ 130 Ce, corroborate this fact. For these reasons, we think that it is needless to account for 
the axial asymmetry in a ’’pure” static study of the equilibrium deformation. However, we have to consider prolate 
(7 = 0°) as well as oblate (7 = 60°) nuclear shapes. In this regard, it is worth remembering that oblate shape given 
by {(3 > 0 ,7 = 60°) is equivalent to the set (/3 < 0 ,7 = 0°). 


A. Comparison between the different contributions entering in the potential energy surface 


It could be useful to compare the importance of the different terms entering in the right hand side of Eq.flG]). 
In this respect, we have drawn in FigH]for axially prolate shape, the four microscopic contributions SE s h e ii(N, (3), 
SEshell(Z, (3), 5P P airing(N, (3), SP pa iring(Z, (3) for the case of 160 Ce as functions of (3. Following the cited order, we 
can say that the difference between the highest and lowest values (in the interval (3 £ [0.0, 0.7]) are respectively about 
11.0MeV,10.5MeV,5.7MeV,3.5MeV for the four corrections. Thus, these variations show that the shell corrections 
5E s heii(N, (3), 5Eshell{Z, /3) are more important than 5P pa iring(N, (3), 5P pa iring(Z, (3) and have a clear minimum at 
respectively f3 = 0.35 and (3 = 0.30. It is well known that for each kind of nucleon the shell correction is in opposite 
phase with respect to the pairing correction (this means for that when 5E s h e ii(N, f3) increases with f3, 6P pa iring (AT, f3) 
decreases and vice versa). 

Contrarily to these curves, the liquid drop model is strictly increasing with f3, and its minimum occurs always at 
the beginning f3 = 0.(spherical shape). When all the contributions are added, the minimum of the potential energy 
surface of the nucleus is reached at about (3 = 0.3 and is mainly due to the shell corrections. When [3 becomes more 
and more, larger the contribution of the liquid drop energy becomes preponderant so that the equilibrium deformation 
occurs generally between (3 = 0 and (3 = 0.4. Because of the convention of the sign stated before, SB m i cro defined 
in Eq. © must be negative in order to increase the binding energy of the nucleus. Since the shell corrections (for 
protons and neutrons) play a major role in 8B m i cro , it is naturally expected that negative (but absolute large) values 
of shell correction contribute to increase the binding energy of the nucleus. In this respect, it is well known that the 
shell correction is essentially determined by the distribution of single-particle levels in the vicinity of the sharp Fermi 
level (defined here as the midway between the last occupied level and the first empty level). Following Ref.[3l|, we 
can state that ’’the nuclear ground state, as well as any other relatively stable state, should correspond to the lowest 
possible degeneracy, or, in other words, the lowest density of state near the Fermi level”. This is illustrated in Fig. 
[2] where the single-particle levels are drawn as function of the deformation (3 (7 being fixed at 7 = 0°). To this end 
we have used the FORTRAN code of Ref.jHJ and [24j. The area where the single-particle level density is low near 
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the Fermi level (black stars) is indicated by a circle. Thus, it is not so surprising that, it is in this region where the 
neutron shell correction becomes the most important, involving a minimum in the PES of the nucleus. 


B. Equilibrium deformations 

They are given in table dm for prolate as well as oblate shapes (see table legend for details). The minima of PES 
for the corresponding wells are denoted minpro and minobl. The deformation energy is defined as the difference 
Edef = Epes{ 0) — Epg S (Jf), i.e., the difference between the potential energy surface for a spherical shape and the 
one corresponding to the absolute minimum of PES. Permanent deformations will be in principle characterized by 
large values of E^ef and are responsible of rotational spectra. 

From this table, some remarks may be drawn: 

(i) Two regions of prolate deformation are found. They occur around TV = 64 and TV = 102 with maximum deformation 
about /3 ss 0.30. The deformation energy (between spherical and deformed shape) is about 6.70 MeV for TV = 64 and 
9.30 MeV for TV = 102 and decreases from either side from these two nuclei. 

(ii) Spherical deformation occur at and near the (magic) numbers TV = 82 and TV = 128 (not shown). 

(iii) The deformation energy decreases from TV = 64 (maximum) to TV = 82 (minimum) and reincreases again to 
TV = 102 (maximum). We have found graphically that the first inflexion point occurs between TV = 72 and TV = 74 
and a second inflexion point is found between TV = 90 and TV = 92. One can consider (somewhat arbitrarily) that 
spherical shapes occur approximately between these two limits. 

(iv) The minima of prolate equilibrium deformations are, by far, always deeper compared to the ones of the oblate 
minima ( minpro minobl). In other words cerium isotopes prefer, by far, prolate shapes. In other words, the 
deformation energy increases in average with the asymmetry 7 . This justifies a posteriori that, in a static study of the 
equilibrium deformation, it is needless to account for axial asymmetry. It is worth to remember that most of nuclei 
of the chart have prolate shape (see Ref. (25j). 

(iv) Even though the experimental deformations are known only in absolute value from B{E 2), a good agreement is 
obtained if one excepts the three ’’nearly magic” nuclei 138 ~ 142 Ce. 

In Fig.® are displayed the present equilibrium deformations, experimental values [l] > our ’’old” calculations Q 
and other studies performed by different authors which are: Kern et al.JHJ, Hilaire and Girod [35), Gotz et a 1. [ 33 J and 
Nix et al. [34j]- All calculations are based on Macro-Micro method (with different mean fields or different parameters). 
Except the one of Ref. [HI which uses Hartree-Fock-Bogoliubov model with Gogny force. 

(j) Near magic number (TV = 82) all calculations give spherical equilibrium deformation whereas experimental results 
are always slightly deformed (even for TV = 82). It seems difficult to overcome this defect with a pure static approach 
which neglects the role of the mass parameters. 

(jj) The overall tendency of these calculations is the same except the fact that HFB calculations differ significantly 
from the others with higher values in some regions. 

(jjj) Apart from HFB calculations, theoretical values are generally quite close from each others. 

(jv) Our old and new calculations give very close results (see Table HIT1) . Thus, even if it is better to choose a proper set 
of mean-field parameters for each nucleus, we do not commit a significant error by taking the same set of parameters 
for nuclei that do not differ strongly by the number of neutrons (TV). 


C. Mass Excesses 

We list from a FORTRAN file (see Fig. [4] ) the results of our theoretical calculations of the binding energies and 
mass excesses (m-excess) for the even-even cerium isotopic chain. For the sake of completeness, experimental mass 
excesses and the ones of the FRDM model (see Ref. 0) are also given. We must keep in mind that only 6 parameters 
are used in the liquid drop model whereas 16 parameters are necessary in the FRDM model. This explains the ’’better 
quality” of the FRDM model. However, we have checked that the variations of binding energy or mass excesses from 
one isotope to the nearest is practically the same in our model and the one of FRDM (the deviations are about 
±0.35MeF). For this reason, the calculation of the two neutron separation energies (see the following subsection 
IIIIDI) will almost be probably the same for the two approaches even though our model is not so accurate. 


D. Transitional regions in cerium isotopes 

In Fig. [5] is shown the gradual transition in the potential energy surface from spherical vibrator to the axially 
deformed rotor when the number of neutrons (TV) increases from 76 to 92. One signature of A'(5) symmetry which is 
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a critical-point of phase/shape transitions (quantum phase transition between spherical and axial symmetries) should 
be a long flatness of the potential energy surface with eventually a weak barrier from prolate to oblate shapes. In 
this figure, for TV > 82, the width of the flatness increases as one moves away from TV = 82 but at the same time 
the difference between oblate and prolate minima and barrier between oblate and prolate shapes also increase. For 
example the differences between oblate and prolate energy minima and barriers for isotopes with TV = 88,90,92 are 
respectively about 1.5MeV, 2.5 MeV and 3.3MeV with energy barrier about 2MeV , AMeV and 5.5MeV respectively. 
The wideness of the bottom of the well must be relativized with the height of the barrier. Thus for the case of TV = 92 
the width is important, i.e. about A/3 « (3 pro — f3 0 u ~ 0.26 — (—0.20) « 0.46 but the barrier is about 5.5MeV and 
therefore seems too high. The case TV = 90 gives a width of A/3 « 0.3 with a barrier of about 4 MeV. For TV < 82, 
the case TV = 76 seems to be relatively equivalent to TV = 90 with a slightly smaller width and a lower height barrier. 
Thus it is difficult to determine clearly the existence of a X(5) critical-point. Thus, everything seems to indicate a 
continuous transition. 

In Fig. [6] is displayed the two-neutron separation energy (TSN) as function of the neutron number TV. A clear 
jump is seen from TV = 82 to TV = 84, i.e. from one major shell to the following. Just before TV = 82 and just after 
TV = 84 the TSN varies more slowly. Far for the ’’jump” the curve becomes quasi-linear. Once again, no special 
behavior is noted around TV = 90 which from Ref.[36| and [37| should constitute with Z « 62 the first order shape 
transition (X(5) critical-point) in the rare earth region. In Ref.[38|, it has been pointed out that ’’Empirical evidence 
of transitional symmetry at the X(5) critical-point has been observed in 150 TVd, 152 Sm, 154 Gd, and 156 Dy”. 

One of the most important signatures of the phase transition is given by a sudden jump in the value of the energy 
ratio R 4/2 = 4^ /2f from one nucleus to the next. We found it useful to compare the experimental values of this ratio 
(see Fig. 0 in the cases of the isotopic chains of Ce and Sm (The experimental values of the considered levels have 
been deduced from the adopted level of ENSDF site [39|). The figure shows clearly two facts. First, the important 
variation of R 4/2 near of the magic number TV = 82 for the both isotopic chains and then, the important difference 
between the behavior the two isotopic chain from TV = 88 to TV = 90. In effect in the case of the Samarium, there 
is a sudden increase of this ratio whereas this is not the case for the Cerium isotopes. This has been attributed to 
the X (5) critical-point symmetry of the nucleus 152 Sm. Thus the present study confirms that cerium isotopic chain 
is characterized by a continuous shape/phase transition. 


IV. CONCLUSION 

Potential energy surfaces have been drawn for the cerium isotopic chain. All even-even nuclei between the two 
drip lines have been considered. To this end, we have used the microscopic macroscopic method in which the quan¬ 
tum corrections have been evaluated by a semi-classical procedure. The microscopic model is based on a ’’realistic” 
Schrodinger equation including a mean field of a Woods-Saxon type. The macroscopic part of the energy is evaluated 
from the liquid drop model using the version of Myers and Swiatecki. The following points must be remembered: 

(i) All equilibrium deformations have been found prolate with an important deformation energy compared to oblate 
shapes. 

(ii) The maximum deformations are of order /3 « 0.3 and are located around TV = 64 and TV = 102 with deformation 
energy about 6 MeV and 9 MeV respectively. The equilibrium deformations decrease as one moves away from these 
two nuclei. 

(iii) Spherical shapes are found in the neighborhood of TV = 82. 

(iv) Good agreement is obtained between theoretical and experimental values if one excepts the area of the shell 
closure TV = 82 where the latter are slightly larger. 

(v) This isotopic chain possesses a continuous shape/phase transition from spherical shapes toward the axially sym¬ 
metric ones. 
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Appendix A: Constants of the binding energy of the liquid drop model 

The constants of Eq. m are defined as follows: 


Cy — ay [l — k/ 2 ] (in the volume term) 
C s = a s [l - nl 2 ] ( in the surface term) 


/ = —-— (relative neutron excess) 

N + Z K 1 

£ = +1 (even — even), 0 (odd), —1 (odd — odd) (in the pairing term) 



5 r 0 




The last correction to the Coulomb energy takes into account that the liquid drop has not a sharp but a diffuse surface 
of the Woods-Saxon type [ ]. The diffuseness parameter is ao and the charge radius ’’contains” ro ( R c h = roA 1 / 3 ). 


Appendix B: Constants of the Woods-Saxon mean potential 


Universal parameters” of the Woods-Saxon central and Spin-orbit potentials entering in Eq[5] 



Central mean field Depth (MeV) 

Central mean field Radius (/to) 

Spin-orbit mean field SO-coupling strength A 35.0 


Neutrons Protons 

Voneut = 49.6(1 - 0.86 1) V 0prot = 49.6(1 + 0.86J) 
Rvneut = 1.347 A 1 / 3 R V prot = 1.275A 1 / 3 


Spin-orbit mean field Radius (/to) 
Central mean Held diffuseness (/to) 
Spin-orbit mean field diffuseness (/to) 


Rso-neut = 1.310A 1 / 3 
a 0 = 0.70 
ao = 0.70 


36.0 

Rso —nprot 

a 0 = 0.70 
a 0 = 0.70 


= 1.200A 1 / 3 


Figures 


Contributions to the Potential Energy Surface of 16 ° 58 Ce 102 



Deformation p 


FIG. 1: Contributions of the shell and pairing corrections for the two kind of nucleons and the one of the liquid drop 

model to the total potential energy surface of the nucleus 160 Ce. 
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FIG. 2: Single-particle energies of the microscopic model as function of deformation for prolate shapes (j3 > 0) for 
the nucleus 160 Ce. Spherical spectroscopic notation is given for spherical deformation (f3 = 0) .The circle in dotted 

line indicates the area of lowest level density. 
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FIG. 3: Theoretical equilibrium deformations for even-even cerium isotopes evaluated by different or similar 

approaches. 
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FIG. 4: Theoretical binding energies and mass excesses of the present approach compared to the experimental mass 
excesses and the ones given by the FRDM model of Ref. []j§]. All energies are expressed in MeV. The experimental 
data as well as the frdrn results have been entered manually in the code. Asterics mean that no experimental data is 

available for the corresponding nucleus. 
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FIG. 5: Shape evolution for cerium isotopes from N = 78 to N = 92. 
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Two-Neutron Separation Energies VS Neutron Number 
For Cerium Isotopes (Z=58) 
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FIG. 6: Two-neutron separation energies (Sin) along the cerium isotopic chain. This quantity is defined as 
S 2 n(A , Z, N) = Bind(A , Z, N ) — Bind(A — 2, Z, N — 2) where the binding energy Bind(A, Z, N) is given by Eq. (JTJ. 
Note that in our approache the neutron drip line (where S-in ~ 0) can be extrapolated around N = 128 for Cerium 

isotopes. 



FIG. 7: I? 4/2 enery ratio as function of neutron number for Cerium and Samarium isotopes. Sudden variations are 
associated with magic closure shells for the both chains ( at IV = 82) and with A' (5) critical point which occurs only 

for Sm (at N = 90). 
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Tables 



CLy 

as 

Cc 

K, C d CLpair 

Myers and Swiatecki 

15.67 MeV 

18.56 MeV 

o 

to 

s; 

cb 

1.79 1.21 MeV 11 MeV 


TABLE I: Parameters of the liquid drop model in the Myers and Swiatecki version [28s] 


N 

A 

/3 P ro minpro /3 0 bi 
(MeV) 

minobl Edef 
(MeV) (MeV) 

l^exp| 

N 

A 

/3 P ro minpro /3 ob i 
(MeV) 

ninobl Edef |/5exp| 

(MeV) (MeV) 

58 

116 

0.30 

0.90 

-0.21 

3.62 

4.80 


92 

150 

0.25 

1.23 

-0.17 

4.45 

5.12 0.31 

60 

118 

0.32 

0.88 

-0.23 

4.07 

5.87 


94 

152 

0.27 

1.21 

-0.19 

5.05 

6.40 

62 

120 

0.32 

1.03 

-0.23 

4.33 

6.19 


96 

154 

0.28 

0.64 

-0.21 

4.94 

7.47 

64 

122 

0.31 

1.16 

-0.23 

4.23 

6.68 


98 

156 

0.29 

0.66 

-0.22 

5.13 

8.44 

66 

124 

0.30 

1.47 

-0.21 

4.15 

6.17 

0.30 

100 

158 

0.29 

0.71 

-0.22 

5.14 

9.08 

68 

126 

0.29 

1.75 

-0.21 

3.87 

5.43 

0.33 

102 

160 

0.30 

0.32 

-0.22 

4.52 

9.27 

70 

128 

0.27 

1.82 

-0.21 

3.48 

4.67 

0.29 

104 

162 

0.29 

0.71 

-0.22 

4.42 

9.08 

72 

130 

0.25 

2.02 

-0.2 

3.27 

3.34 

0.26 

106 

164 

0.29 

1.00 

-0.23 

4.23 

8.44 

74 

132 

0.20 

1.90 

-0.17 

2.60 

1.97 

0.26 

108 

166 

0.28 

1.16 

-0.23 

3.92 

7.57 

76 

134 

0.16 

1.28 

-0.14 

1.63 

0.93 

0.19 

110 

168 

0.27 

1.46 

-0.21 

3.84 

6.39 

78 

136 

0.10 

0.04 

-0.07 

0.18 

0.19 

0.17 

112 

170 

0.25 

1.68 

-0.20 

3.55 

5.33 

80 

138 

0.00 

-1.93 

0.00 

-1.93 

0.00 

0.13 

114 

172 

0.25 

1.97 

-0.19 

3.20 

4.19 

82 

140 

0.00 

-3.96 

0.00 

-3.96 

0.00 

0.10 

116 

174 

0.2 

1.93 

-0.17 

2.79 

2.95 

84 

142 

0.00 

-2.07 

0.00 

-2.07 

0.00 

0.12 

118 

176 

0.17 

1.71 

-0.16 

2.17 

1.68 

86 

144 

0.15 

0.02 

-0.06 

0.53 

0.50 

0.17 

120 

178 

0.14 

1.39 

-0.14 

1.60 

0.55 

88 

146 

0.19 

0.73 

-0.11 

2.43 

1.99 

0.17 

122 

180 

0.0 

0.3 

0.00 

0.30 

-0.15 

90 

148 

0.23 

1.15 

-0.14 

3.76 

3.15 

0.25 

124 

182 

0.0 

-1.08 

0.00 

-1.08 

-0.08 


TABLE II: Equilibrium deformations as well as deformation energies for the cerium isotopic chain. The columns 
give successively the number of neutrons (N), the mass number (A), the prolate equilibrium deformation (/3 pro ), the 
minimum of the prolate well ( minpro ), the oblate equilibrium deformation (/ 3 0 u ), the minimum of the oblate well 
( minobl ), the deformation energy (Edef ,see text), the experimental equilibrium deformation (/3 eX p)- Note: The 
deformation energy is always given for the prolate equilibrium shape because no absolute minimum is obtained for 

oblate shape. 


Cerium (Z = 58) 

N = 58 

60 

62 

64 

66 

68 

70 

72 

74 

76 

78 

80 

82 

Present (3 

+0.30 

+0.32 

+0.32 

+0.31 

+0.30 

+0.29 

+0.27 

+0.25 

+0.20 

+0.16 

+0.10 

+0.00 

+0.00 

Old /3 

+0.28 

+0.30 

+0.31 

+0.31 

+0.31 

+0.30 

+0.27 

+0.24 

+0.22 

+0.18 

+0.06 

+0.11 

+0.00 

Present Edef (MeV) 

4.80 

5.87 

6.19 

6.68 

6.17 

5.43 

4.67 

3.34 

1.97 

0.93 

0.19 

0.00 

0.00 

Old E def (MeV) 

4.82 

5.77 

6.03 

6.31 

7.08 

5.36 

4.41 

3.35 

2.13 

0.77 

0.00 

0.24 

0.00 


TABLE III: New equilibrium deformations and deformations energies vs old @. 
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